The Column Density Distribution of the Lyman-Alpha Forest: 
Measure of Small Scale Power 
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Absorption lines in the Lya forest can be thought of as peaks in neutral hydrogen 
density along lines of sight. The column density distribution (the number density of 
absorption lines as a function of column density) is then a statistic of density peaks, 
which contains information about the underlying power spectrum. In particular, 
we show that the slope of the distribution provides a measure of power on scales 
smaller than those probed by studies of present-day large scale structure. 



Two examples of power spectra are shown in Fig. la. Plotted are cro(ks) = 

*£_ 

[/ 4:irP(k)e k s 2 k 2 dk] - 5 versus ks, where oo(fcs) can be regarded as the 
amount of power on scale ks for the power spectrum P(k), linearly extrap- 
olated to the redshift of interest. The CHDM model, because of neutrino free 
streaming, has less power than the CDM model on the scales shown. What 
kind of column density distribution would each predict? 

Several steps are involved in the computation. We briefly mention a few 
important ones and refer the reader to Hui, Gnedin & ZhangEI for details. 

Let us denote the local gas overdensity by 8b- The neutral hydrogen density 
riHi is then proportional to (1 +6b) 2 T~°''J^ l assuming ionization equilibrium, 
where T is the temperature and Jhi is the intensity of the radiation background 
at the hydrogen ionizing frequency. Using the approximate relation T — Tq(1 + 
8b) a where To is the temperature at mean gas density and a is some prescribed 
power law index (see Hui & Gnedin 0), one deduces that nni is proportional 
tofigT -°- r j£(l + *6) a -°- 7a - 

Now, consider a density peak along a line of sight, its neutral hydrogen 
column density is given by: 
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where r is the proper distance along the line of sight. We have Taylor expanded 
the density field around the peak(i.e. rewriting 1 + 8b as e^ and expanding £ 
to second order) so that the integrand uhi becomes a Gaussian. 

Two important conclusions follow from this simple expression. First, for 
given 8b and its second derivative, changing the combination of cosmological 
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parameters Q%Tq J^ , which is uncertain by a factor of about 10 observa- 
tionally, has the effect of rescaling JVhi • This means that for a column density 
distribution approximated by a power law (see Fig. lb), the 'left-right' place- 
ment of the distribution is uncertain by an order of magnitude; in other words, 
at a fixed column density, this means a rescaling in the normalization of the 
distribution. Because of this it is hard to use the normalization of the column 
density distribution to learn something about the power spectrum. 

Second, the column density depends on 1 + <5& and its second derivative. 
This means that if one knows the probability distribution of 1 + 8b and its first 
and second derivatives (with the first derivative vanishing) at any given point, 
one can calculate the column density distribution. To do so, one would need a 
way of evolving the density field 8b. We use the Zel'dovich approximation as an 
efficient method which is sufficiently accurate in the mildly nonlinear regime 
8 b ^ 5, corresponding to TVjji ~ 10 145 cm 2 at z = 3 for typical cosmological pa- 
rameters. An appropriate smoothing scale has to be chosen when using the ap- 
proximation, to take into account orbit-crossing and also smoothing due to gas 
pressure on small scales. Without going into details^it turns out for most inter- 
esting cosmological models, the scale lies in the range IMpc -1 <fc<10Mpc _1 
where k is the comoving wave number. (Note this is where Co ~ 1 in Fig. 
la). This is relatively small scale compared to scales familiar in studies of 
present-day large scale structure (0.01 Mpc -1 <k <1 Mpc -1 ). 

Results of computations using the Zel'dovich approximation are shown in 
Fig. lb. A comparison is made with the results from a full hydrodynamic 
simulation and the agreement is encouraging. We note that the CHDM model 
produces a steeper distribution compared to the CDM model. This can be 
understood qualitatively by recalling that the CHDM model has less power 
than the CDM model on the relevant scales (Fig. la). This means that by 
z = 3, the CHDM density field is relatively less nonlinear compared with 
the CDM one i.e. there are proportionally fewer high density peaks 
compared to intermediate density ones, hence the steeper column 
density distribution of the CHDM model, fl -fl 



d There are also other subtleties such as the effects ot-peculiar velocities, the definition of 
peaks, etc for which the reader is referred to Hui, et al.tj 

e Note also that, if one goes to sufficiently low density, because the CHDM has not de- 
veloped very underdense regions, there would also be proportionally fewer very low density 
peaks compared to intermediate ones, but the column density range we show is not low 
enough to unambiguously see this effect. 

' The above argument assumes that high column density also means high density 1 + <5(, 
(or high peak height) and vice versa. From Eq. ll it is clear column density also depends 
on the peak width (the term involving second derivative of density). However, peak width 
and height are correlated and the correlation works out so that column density still increases 
with peak height. 
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Figure 1: la on the left shows cro versus ks as defined in the text, for two power spectra 
at z = 3: solid line is a CDM (Cold Dark Matter) model with erg = 0.7 and dashed line is a 
COBE-normalized CHDM (Cold+Hot Dark Matter) model with U„ = 0.2, a tilt of n = 0.9 
and no tensor modes. Both have Hq = 50kms — 1 Mpc — 1 and fi = 1. lb on the right shows 
the resulting column density distributions (number of absorption lines Nn nB per unit column 
density JVhi per unit redshift z versus column density), open triangles for CDM and open 
squares for CHDM, computed using the Zel'dovich approximation. Solid triangles represent 
the same CDM model but computed using a full hydrodynamic simulation. Two different 
values of Sl^T^ ' J^ are chosen for the CDM and CHDM models respectively, which shifts 
the normalization of the points but not the slope, to clearly show the difference in their 
predicted slopes. See Hui, et al. x for details of the computations. 



To conclude, while the normalization of the column density distribution 
provides only a weak constraint on the power spectrum due to large observa- 
tional uncertainties in fi^TJ" ' J^ T , its slope can provide useful information 
on the amount of power on scales complimentary to those of other traditional 
cosmological studies. 

This work was supported in part by the DOE and the NASA under grant 
NAG 5-2788. 



References 



1. L. Hui, N. Y. Gnedin, Y. Zhang, preprint astro-ph 9608157 

2. L. Hui, N. Y. Gnedin, preprint |astro-ph 9612232| 



